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■ Abstract 

r~| , We show that the g-Digamma function "0<j f° r < q < 1 appears in 

an iteration studied by Berg and Duran. In addition we determine the 
probability measure v q with moments 1/Ylk=i(l ~ <?)/(! ~~ Q k )i which are 
g-analogues of the reciprocals of the harmonic numbers. 
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1 Introduction 

For a measure fi on the unit interval [0, 1] we consider its Bernstein transform 



X 



b(a*)W = jT \^Mt), (i) 

as well as its Mellin transform 

M(fi){z) = [ t z d^{t), Uz > 0. (2) 
Jo 

These functions are clearly holomorphic in the right half-plane ^tz > 0. 

The two integral transformations are combined in the following theorem from 
[3] about Hausdorff moment sequences, i.e., sequences (a„)„>o of the form 

a n = [ t n dfi(t), (3) 
Jo 

for a positive measure n on the unit interval. 
* Corresponding author 
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Theorem 1.1 Let (a n ) n >o be a H aus dor ff moment sequence as in ([3]) with \i ^ 0. 
Then the sequence (T(a n )) n > defined by T(a n ) n = l/(a + • • • + a n ) is again a 
Hausdorff moment sequence, and its associated measure T(fx) has the properties 
f(/i)({0}) = and 

B{fi){z + l)M(T(n))(z) = 1 for 3fc > 0. (4) 

This means that the measure T(^) is determined as the inverse Mellin transform 
of the function 1/B(fi)(z + 1). 

It follows by Theorem 11.11 that T maps the set of normalized Hausdorff mo- 
ment sequences (i.e., a$ = 1) into itself. By Tychonoff's extension of Brouwer's 
fixed point theorem, T has a fixed point (m n ). Furthermore, it is clear that a 
fixed point (m n ) is uniquely determined by the equations 

(1 + mi + . . . + m n )m n = 1, n > 1. (5) 

Therefore 

-^1 + ^-1 = 0, (6) 

m n 

giving 

-1 + V5 - V5 - 1 
mi = , m 2 = 

Similarly, T maps the set M*([0, 1]) of probability measures on [0, 1] into 
itself. It has a uniquely determined fixed point u and 

m„= f t n du(t), n = 0,l,.... (7) 
Jo 

Berg and Duran studied this fixed point in |1],[5], and it was proved that 
the Bernstein transform / = B(u) is meromorphic in the whole complex plane 
and characterized by a functional equation and a log-convexity property in ana- 
logy with Bohr-Mollerup's characterization of the Gamma function. Let us also 
mention that u has an increasing and convex density with respect to Lebesgue 
measure m on the unit interval. 

An important step in the proof is to establish that u is an attractive fixed 
point so that in particular the iterates T on (Si) converge weakly to u>. Here and 
in the following 5 a denotes the Dirac measure with mass 1 concentrated in a G 1. 

It is easy to see that T(5i) = m, because 

T(l,l,...) n =_j—= [ t n dt. 
n+1 Jo 

It is well-known that the Bernstein transform of Lebesgue measure m on 
[0, 1] is related to the Digamma function if), i.e., the logarithmic derivative of the 
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Gamma function, since 



^^ rft = ^ +1 ) +7 = g__, S2>0 , (8) 



cf. 8.36]. Here 7 = — is Euler's constant. 
Therefore v\ := T(m) = T o2 (Si) is determined by 

M{vi){z) 



B(m)(z + 1) tjj(z + 2) +7' 
The measure V\ = T(m) is given explicitly in [3] as 



ui = (J2<xnt in )dt, (9) 



n=0 



where £q = 0, G (n, n + 1), n = 1, 2, . . . is the solution to ip{l — = —7 and 
ct„ = — £ n ). The moments of the measure v\ are the reciprocals of the 

harmonic numbers, i.e., 

1 1 /n+l-.X- 1 

The purpose of this paper is to study the first elements of the sequence T on (5 q ), 
where < q < 1 is fixed. The reason for excluding q = is that T(5q) = S\. 
Since u is an attractive fixed point, we know that the sequence converges weakly 
to OJ. 

The first step in the iteration is easy: 

00 

= (l-g)£gV (11) 

because 

/ t* d?(5 q )(t) = ^r = (i-g)E^ z . (12) 

This shows that T(5 q ) is the Jackson d q t-mea.sure on [0, 1] used in the theory 
of q- integrals, cf. [8]. It is a g-analogue of Lebesgue measure in the sense that 
d q t — > m weakly for q — > 1. 

It is therefore to be expected that v q := T(d q t) = T o2 (6 q ) is a g-analogue of 
the measure Ui, and we are going to determine v q as closely as possible. We have 

MM{Z) = jlkry (13) 
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where f q is defined as the Bernstein transform of d q t: 

— d «t = a - 9) [z + E 9* J • ( 14 ) 

This formula is a g-analogue of (jHJ). The moments of v q are g-analogues of (|T0l) 

Our main result is the following (note that the Haar measure on the multi- 
plicative group ]0, oo[ is dt/t): 

Theorem 1.2 The measure v q has a continuous density u q {t) with respect to 
dt/t on ]0, 1]. It is C°° on each of the open intervals ]q +1 , q k [, k = 0, 1, . . . with 
jump of the derivative of size q k / (1 — q k )(l — q) at the point q h ,k = 1,2,.... 
Furthermore, lim^o v q {€) = 0. 

Remark 1.3 It follows that the behaviour of u q (t) is oscillatory, and therefore 
quite different from that of lim^i u q (t), which is increasing and convex. In fact, 
it follows from (Q that 

oo 

hmvJt) = Va/ +1 , < t < 1. 

n=0 

See Figure 1 and 2 which shows the graph of v q (e~ l ) for q = 0.5 and q = 0.9. 



2 Proofs 

Jackson's g-analogue of the Gamma function is defined as 

T q {z) = ^^{l-qf-\ 

cf. [8], and its logarithmic derivative 

d °° q k+z 

ip q {z) = -j- log T q (z) = - log(l - q) + log q 7 ~ fc+z ( 16 ) 

Z k=0 y 

has been proposed in [TT] as a g-analogue of the Digamma function ip. See also 
the recent paper [12]. We define the g-analogue of Euler's constant as 

oo u 

7 9 = -ip q (l) = log(l - q) - logg^ fc . (17) 

fc=i ^ 
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The Bernstein transform f q of d q t is given in (jHJ), hence 



n / \ OO OO 

" k=\ n=0 

oo / oo 

fc(n+l) _ fc(n+l+z)^ 



+ E 5> fe(n+1) - ? A 



= z 

n=0 \k=l 

which shows the close relationship with the g-Digamma function. We will be 
using another expression for f q (z)/(l — q) derived from (TH|) . namely 



/«(*) . V- q k 



Eq 

with 

00 q' 

q 



kz 



q 

\2 



'-Ei ( 19 ) 



Clearly, g/(l — g) < c q < g/(l — q) 2 for < g < 1 and g H- c g is a stricly increasing 
map of ]0, 1[ onto ]0, oof. We mention two other expressions 



00 

n=l n=l 



J>(n)g" = ^(l-(g";g) c 



where d(n) is the number of divisors in n, see [7J p. 14]. 

In order to replace the Mellin transformation by the Laplace transformation 
we introduce the probability measure r q on [0, oof which has v q as image measure 
under t —> e _t , hence 



00 



C{r q ){z)= I e- tz dr q (t) 1 



/«(*+!)" 

The analogue of Theorem 11.21 about the measure r q is given in the next theorem, 
which we shall prove first. 

Theorem 2.1 The measure r q has a continuous density also denoted r q with 
respect to Lebesgue measure on [0,oo[. It is C°° in each of the open intervals 
]nlog(l/g), (n + 1) log(l/g)[, n — 0, 1, . . . with jump of the derivative of size 

In 

at the point n log(l/g), n = 1,2,.... Furthermore, lim^oo T q (t) = 0. 
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Proof of Theorem \2.1\ Introducing the discrete measure 



00 2k 



k=i 



of finite total mass 



we can write 



fAh = c i -v/i 1 -?) < c g> ( 21 ) 

1 + Cg + Z ~ C(I*)(Z), 



1-q 



hence 



Let p g denote the following exponential density restricted to the positive half-line 

p q (t) = exp(-(l + c q )t)Y(t), 

where Y is the usual Heaviside function equal to 1 for t > and equal to zero 
for t < 0. Its Laplace transform is given as 

e- tz p q (t) dt = (1 + c q + z)-\ 

but this shows that (}2"2"]) is equivalent to the following convolution equation 

00 00 
(1 - q)r q = p q *J2(»*P< i y n = Yl Pl (n+1) * V* n - (23) 

n=0 n=0 

This equation expresses a factorization of (1 — q)r q as the convolution of the 
exponential density p q and an elementary kernel iV* n with N = p * p q . For 
information about the basic notion of elementary kernels in potential theory, see 
[6j p. 100]. All three measures in question r q ,p q and J2o^(p * Pg)* n are potential 
kernels on K in the sense of [6]. 

The measure p* n ,n > 1 is a discrete measure concentrated in the points 
k log(l/g), k — n, n + 1, . . .. The convolution powers of p q are Gamma densities 

-in 

as is easily seen by Laplace transformation. 

Clearly, p q * p is a bounded integrable function with 



\\Pq*V\\oc < ||Pff||oo|H|l < Cq, \\Pq*P\\l=\\Pq\\l\\l*\\l<T-T—, ( 24 ) 

-L ~r C q 
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and then p q * (p q * p)* n , n > 1 is a continuous integrable function on M, vanishing 
for t < n\og(l/q) and for t — > oo. Furthermore, 

\\p q * (p q * p)* n \\oo < (C q /(1 + C q )) n , 

and this shows that the right-hand side of (1231) converges uniformly on [0, oof, so 
(1 — q)r q has a continuous density on [0, oof tending to at infinity. 
For n > 1 and x G [n log(l/g), oof we get 

p* q {n+1) * p* n {x) 

x (x - t) n e _ {1+Cq)(x _ t)Y _ 
n\ 

= e -(i+c,)« y (x-fclog(l/g))^ _ fc(1+Cg)y -k\og{l/q))p? n {k\og(l/q)) 

which is a finite sum, and 

/^(Mog(l/g)) = £ JJ^—-, Ar = n,n+1,.... 

Pi+...+p„=fc j=i 



q 2 x 



In particular, 

/i*>log(l/g)) , 

1 - g 

For n > and < x < (n + 1) log(l/g) we then get 
(1 - q)r q (x) = (25) 

e -(i + c)* q'^Yix - j log(l/g)) {x - jl °f 1/q))k ^(j Iog(Vff))- 
i=o fc=o 

On [0, log(l/g)[ it is equal to exp(— (1 + c g )x), on [log(l/g), 2 log(l/g)[ it is equal 

to 

exp(-(l + c q )x) + f^O* - log(l/g))J , 
on [2 log(l/g), 31og(l/g)[ it is equal to 

exp(-(l + c q )x) (1 + y^(x - log(l/g))+ 



Q 



2(l-c g ) 



x - 2 log(l/g)) + -(x - 2 log(l/g)) 2 , 



1-g 2 v m 2(1 -q) 2 

etc. 

Using the expression (|2"5"1) it is possible to calculate the derivative of (1 — q)r q 
from the right and from the left at the point x = n\og(l/q),n > 1. The difference 
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between the right and the left derivative equals q 2n /(l — q n ) and this gives the 
jump J„ of pop . □ 

It is straightforward to transfer the results of Theorem 12.11 to give Theo- 
rem 11.21 using that v q is the image measure of r q under t i— > e _i , hence u q (t) = 
r 9 (log(l/t)), 0<t<l. 

Remark 2.2 The representation f )22|) and Theorem 12. II show that (1 — q)r q is a 
standard p-function in the terminology from the theory of regenerative phenom- 
ena, cf. [ID] . 




Figure 1: The graph of (1 — q)r q on [0, 3 log(l/g)] for g = 0.5 



3 Further properties of r q 

Formally, by Fourier inversion we get that 



2vr 7_oo + iy) 
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Figure 2: The graph of (1 — q)r q on [0, 3 log(l/g)] for q = 0.9 



The function l/f q (l + iy) is a non-integrable L 2 -function, so the formula holds 
in the L 2 -sense. To see this we notice that 

^l = l-q+ [ e - tz h q {t)dt, Uz>0, 
z Jo 



where 



,(*) = (l~q) TT 



(26) 



fc>i/log(l/g) 

In particular 

ML±M = i _ q + r e -^e-%(t) dt, 
1 + W Jo 

and since e~ l h q (t) is integrable, it follows from the Riemann-Lebesgue Lemma 

that we get the asymptotic behaviour 



f q (l + iy) ~ (1 -q)(l + iy), \y\ -too. 



(27) 



Furthermore, we notice that 



Uf q (l +iy) = l- q +(l- q )J2 -JL_ (1 - q k cos(% log(g))) 



k=l 

9 



hence 

oo 

i<K/ q (l + ^)<l-g + ^/(l + g fc ), 



k=l 



showing that 3?/ 9 (l + iy) is bounded below and above. It follows that the sym- 
metrized density 



r q (x) if x > 0, 



<p q {x) = { (28) 
t q (—x) if x < 0, 

is the Fourier transform of the non-negative integrable function 

l/,(i + <v)l 2 ' 

and therefore <p q (x) is continuous and positive definite, so r q is the restriction to 
[0, oof of such a function. 



Remark 3.1 The function f q defined in (114]) is a Bernstein function in the sense 
of [6], but not a complete Bernstein function in the sense of [13], because f q (z)/z 
is not a Stieltjes function as shown by formula ( 126]) . This is in contrast to 

lim/,(2) = ip(z + 1) +7, 
g->l 

which is a complete Bernstein function, cf. [3]. 

4 Relation to other work 

The transformation T can be extended from normalized Hausdorff moment se- 
quences to the set JC = [0, 1] N of sequences (x n ) = (x n ) n >i of numbers from the 
unit interval [0,1]. This was done in [2J, where T : /C — > K. is defined by 

< r <*"»"- 1 + :tl+ 1 ... + x„ - " SL < 29 > 

The connection is that a normalized Hausdorff moment sequence (a«) n >o is con- 
sidered as the element (a n ) n >i G /C. 

Since T is a continuous transformation of the compact convex set /C in the 
space R N of real sequences equipped with the product topology, it has a fixed 
point by Tychonoff's theorem, and this is (m n ) n >i. 

There is no reason a priori that the fixed point (m n ) of ( |29l) should be a 
Hausdorff moment sequence, but as we have seen above, the motivation for the 
study of T comes from the theory of Hausdorff moment sequences. 
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Although T is not a contraction on /C in the natural metric 

oo 

d((a n ),(b n )) = ^2~ n \a n - b n \, (a n ),(b n ) G K, 

n=l 

it was proved in [2] that T maps /C into the compact convex subset 

C = {(O e/C | ai > §}, 

and the restriction of T to C is a contraction. It is therefore possible to infer that 
(m n ) is an attractive fixed point from the fixed point theorem of Banach. 
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